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1 Introduction 

This note is for an introduction to the validated numerics study of values of 
zeta and L - functions. The subject matter of this talk lies in the area between 
Analysis, Algebra and Validated Numerics |l], [7). More specifically, I want to 
discuss the relations between the algebra and analysis of zeta and L - functions 
and the evaluation of values of the functions. The talk is on some observations 
and partial results related to more general cases, and I would like to ask the 
permission of the listeners for replacing the notes of my talk by a very introduc- 
tory exposition. The reason for this "going down" and taking very elementary 
examples as the object for discussion lies on my side, i.e. it was useful for the 
study of more general cases to review the simplest cases. 
Why Zeta and L-functions are considered? 

As noted by Don Zagier, zeta and L-functions of various sorts are all-pervasive 
objects in modern number theory and its applications ||. They appear in num- 
ber theory, in algebraic geometry, in algebraic K-theory, in representation theory 
and in the theory of automorphic forms. Recently zeta and L-functions appear 
in the theory of dynamical systems ||l8| , |l9| . Many mentioned applications con- 
nected with computations of zeta and L-functions values in "critical" points 
(following Delignc [fiof an integer m is called critical if neither m nor h — m is 
a pole of 7(s), where 7(3) is a "gamma factor", which equal to an exponential 
function A s times a finite product of terms T(l/2(s + m)) such that the product 
L*(s) = 7(s)L(s) satisfies L*(s) = wL*(k — s) for some integer h > and sign 
w = ±1). Then the corresponding critical value L(m) should have the form 

L(m) — A(m)Q(m) (m critical) 

where f2(m) is a predictable period (the integral over some algebraic cycle 
of an algebraic differential form defined over a number field) and A(m) is an 
algebraic number belonging to a predictable number field. 
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The organization of this note is as follows: In Section 1 we remind definitions and 
some properties of Riemann zeta function, Dirichlct L- functions, Dedekind zeta 
function, Artin-Hasse zeta functions and Hasse-Weil zeta function of elliptic 
curves. Section 2 comments some applications of zeta and L— functions. In 
Section 3 we consider some validated numerics aspects of zeta and L— functions 
and interval evaluation of Riemann zeta function for Re{s) > 1 and values of 
some Dirichlet L- functions at point 8 = 1. Also we consider some aspects of 
interval evaluation of Hasse-Weil zeta functions of elliptic curves at point s = 1. 



2 Zeta and L - functions 

Let us remind definitions and some properties of Riemann zeta function, Dirich- 
let L-functions, Dedekind zeta function, Artin-Hasse zeta functions and Hasse- 
Weil zeta function of elliptic curves. 



2.1 Riemann's zeta function 

Let 



V — = Y\(l - p~ s r\n G N, p = 2, 3, 5, 7, . . . (primes) (1) 
^— ' n s - L± 

n p 

be the Euler product. Euler used this formula principally as a formal identity 
and principally for integer values of s. Let 



c(*) = E^- ( 2 ) 



Dirichlet, unlike Euler, used the formulae (1.1), (1.2) with s as a real variable, 
and, also unlike Euler, he proved rigorously that (1.1) is true for all real s > 1. 
Riemann defines the function [l2[ 

+oo 

_ rc-fc+i» I JteT ^ (3) 

Im J exp(x) — 1 x 

+oo 

+oo 

where J ... is a counter integral. The limits of integration are intended 

to indicate a path of integration which begins at +oo, moves to the left down 
the positive real axis, circles the origin once in the positive (counterclockwise) 
direction, and returns up the positive real axis to +oo. If ((s) is defined by 
formula (1.3), then for real values of s greater than one, ((s) is equal to function 
(1.2). The formula (1.3) defines a function £(s) which is analytic at all points 
of the complex s— plane except for a simple pole at s = 1. The function £(s) 
satisfies the functional equation 

C(s) =r(-(s + l))(2^) s - 1 2sin(s7r/2)C(l-s). (4) 
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Let Bk be the fc-th Bernoulli number in even numeration Bq — 1,-Bi = — i, . . . . 
For n € Z Euler have found that 

C (2») = ^) 2 "(- 1 )" +1 ^« 



2(2n)! 

If n > then £(— 2ri) = 0. Functional equation (1.4) gives for m E N that 

C(l-2m) = (-l) 1 - 2m |^Ii. 

2.2 Dirichlet L-functions 

Let p be a prime. Let uj 13 ^ 1 = 1. Let g be a primitive root (mod p). Let 
z^(n) be the index of natural n relatively to the primitive root g, i.e. g"W = n 
(mod p). 

Elementary Dirichlet character to a (mod p) is the function 

x(n)=oj u{n \ x(n) = 0ifp| n. 

There are tp(p) = p — 1 distinct functions: x( n ) = e p ~ 1 j 1 < to < p — 1. Let s 
be a positive real. Let 

oo 
n=l 

The function L w (s) is called a Dirichlet L- function. 

Let x( n ) be a Dirichlet character to a (mod q). In the case there exists a 
Dirichlet L-function L{s,x) §]• Let e q {m) = e « . Let 



m— 1 

be a Gauss sum. Let now ^ be a primitive character, and let || 

0, if x (-i) = i 

1, if%(-l) = -l 
Let 

ZM = (-)(- 1 / 2 )^r(i(. s + a))L( s , x ). 

q 2 

Then 

€ (1 X ) (6) 

is the functional equation for L(s,x). Dirichlet L— functions are a partial case 
of Dirichlet series. 
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2.3 Dedekind zeta functions 

Dedekind extended the Riemann zeta function from natural numbers n to all 
integer divisors (ideals) of an algebraic number field if. The Dedekind zeta 
function Ck{s) is defined by the formula 

a v ' 

Here a runs over all integer divisors (ideals) of the field if and N(a) is the norm 
of a. 

Let m € N, £ m = 1. A field if is the m-cyclotomic field if if = Q(£), £ m = 1. 
Let m = p, where p is a prime. In the case 



c 



Here G(s) = 0p| P (^ — TfTpp ) 1 anc ^ P runs over a ^ P r i me divisors of the field 
Q(-v / T), p | P- Let x be all Dirichlet characters (mod p). 
Proposition. |Q Suppose if is the p~ cyclotomic field. Then the Dedekind's 
zeta function of if is the product of G(s) and tp(p) Dirichlet L— functions runs 
over all Dirichlet characters x '■ 

Ck(s)=G(s)1[L(s,x). 



Let p be a positive integer, p = 1 (mod 4), and if = Q( v / p). From Siegel's 
follows, that 

l<b<y/p, bodd 

where a\(n) is the sum of the divisors of n. 

2.4 Zeta and L-functions of elliptic curves 

p3f Let -B/Q be an elliptic curve given in Weierstrass form by an equation 

E : y 2 + a\xy + a^y = x 3 + ci2X 2 + a^x + ae, (8) 

and let 62, 64, be, C4, C6, A, j be the usual associated quantities p3[ . Let 
now (||)be a global minimal Weierstrass equation for £7 over Z. For each prime 
p the reduction (||) (mod p) defines a curve E p over the prime field F p . Let 
A p denote the number of points of E p rational over F p . Let 

t p = 1 +p - A p . 
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If p /fA, then t p is the trace of Frobenius and satisfies \T V \ < In the case 

Artin-Hasse zeta function of the elliptic curve Ep is i 



i - t pP -° + p 1 - 28 
^^ = (i- p -)(i-pi-') - (9) 



If p | A, then £J P is not an elliptic curve and has a singularity S. In the case 

!0, if S is a cusp, 
1, if S is a node, 
— 1, if S is a node with tangent quadratic over F p 

The Hasse-Weil L— function of E/Q is defined by equation 

M s ) = IT n 7 II i ; =tt i-2, ( 10 ) 



3 Evaluation of values of zeta and L - functions 
and their applications 

Many results in number theory are based on evaluation of values of Riemann 
zeta and Dirichlet L-functions (zeros of zeta and L- functions, prime number the- 
orem, distribution of primes in arithmetic progressions, density theorems). Here 
we remind some another applications (please, see also very interesting paper S). 

Evaluation of the index of elliptic operator and zeta functions [^o[ || 

Let A be a self adjoint positive pseudo-differential (Calderon-Zygmund) op- 
erator. Let A has no eigenvalues A on negative real axis. We can define 

A* = ±-J\ s {A-\)- l d\, 

where T is the counter which begin at — oo, moves to the right up the negative 
real axis, circles the origin once in the clockwise direction, and returns down 
the negative real axis to — oo. 

Let A = —-^i + P be the operator on unit circle S 1 , Pf = ^ f(x)dx. 
Eigenvalues of A are natural numbers 1, 1, 1, 4, 4, ... , m 2 , to 2 , . . . and if Xj is the 
j eigenvalue of A, then trA s = Z{s) = 1 + X j = 1 + 2 C(-2s). Suppose 

now that D is an elliptic operator and D* is the adjoint of D. Then 

A = 1 + D*D, Ai = 1 + DD* 



5 



are positive self-adjoint operators. Let 

Q{s) =trA7*,i = 0,l. 
Then the difference of £o and Ci gives the index of the operator D. 

The symplectic volume of moduli space El 

Let A4(n, d) be the moduli space of (semi)stable holomorphic vector bundles 
of rank n, degree d and fixed determinant on a compact Riemann surface S. 
The symplectic volume of Ai(2, 1) is given by 

vol(M(2, !)) = (! 1 ^ (2 ^ 2) 



22g-3 I 2^-2^29-2 ' 

Here g is the genus of £. If n,d,g are sufficiently large, then an evaluation of 
vol(M.(n, d)) is desirable. 

Generalized Riemann conjecture 

Let s = a + it. Let L(s, \) be a Dirichlet L— function. Then 

L( S , X )^0ifa> I (11) 

The Hilbert modular group and the volume of orbifolds [jl6| 

Let K be a, real quadratic field over Q and R the ring of algebraic integers 
in K. Let H be the upper half plane of C. The group SL(2, R) acts on H x H. 
The Hilbert modular group L = SL(2,R)/(1, — 1) acts on H x H effectively. 
The volume of X = T\(H x 77) is given by 

uo/(X) = 2Cx-(-l). 

The Birch and Swinnerton-Dyer conjecture for elliptic curves [E3| 
1,1 

The Hasse-Weil L— function (|l^) converges for Re(s) > I because it is dom- 
inated by the product (£(s — ^)) 2 - 

The Birch and Swinnerton-Dyer conjecture for elliptic curves is 

BSDC1 Let r be the rang of E(Q). Then L(E, s) has a zero of order r at s = 1. 

In particular, /^(-E 1 , s) has the Taylor expansion 

L(E, s) — Ce(s — l) r + higher terms with some constant Ce ^ 0. 

The second Birch and Swinnerton-Dyer conjecture for elliptic curves concerns 
this constant. Arithmetic investigations involve derivatives of L(E, s) and com- 
putation of their values M] . 
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4 Validated numerics of values of zeta and L - 
functions 



4.1 Validated numerics and zeta functions 

Remind at first the interval interpretation of inequality f(z) ^= 0. It appears in 
formula ( pi] ) and in many another cases in the theory of zeta and L— functions. 
The inequality has the interpretation: 

(i) If f(z) € R then (f(z) > 0) V (f(z) < 0) and f(x) € A for a real interval 
A, A. 

(ii) If f(z) 6 C, /(z) ^ R and we use complex rectangular arithmetic 7C 
then (Re(z) > 0) V (iZe(z) < 0) V (Im(z) > 0) V (im(z) < 0). So /(z) G A, A e 
IC, ^ A. Side by side with rectangular arithmetic there are complex circular 
arithmetic H and complex sector arithmetic Jl7| . In some cases they are very 
useful for evaluation of f(z) ^ 0. But as this requires the consideration of nu- 
merous different cases we omit it. In preceding sections we have reminded some 
results about values of zeta functions at integer (critical) points. A common 
problem of validated numerics to find optimal interval evaluation of the value 
of a function with interval variables and parameters. In many cases for the 
purpose we can use: 

(i) programming languages and compilers for scientific computations 
(Pascal-XSC, ACRITH-XSC, Fortran-XSC, Oberon-XSC); 

(ii) computer algebra systems with interval packages; 

(iii) highly accurate, extended computer arithmetics with standard 
and special functions for real and interval data Q . 

H. Stetter gave in |22| and in his another papers applications of computer al- 
gebra in interval methods. In the frame of zeta and L— functions we can apply 
in some cases algebraic-analytic approach and corresponding computer algebra 
for exact presentation of zeta or L— function as a finite sum and a remain- 
ing term. Then the sum and the remaining term are interval evaluated plf . 
Computer algebra systems include now in their standard library functions some 
kinds of Zeta functions. For instance Maple V Release 4 includes JacobiThetal- 
Theta3 (Theta functions are connected with Zeta functions), WeierstrassZeta, 
Riemann Zeta and Hurvitz Zeta. So we can in most interesting cases compute 
(not guaranteed) values of the functions. In some cases the Maple computing 
Zeta(n,a,m),n € Z, a, m £ Q in closed form (by n, a, In). Standard functions 
and packages of Maple give possibility to compute values of some Zeta and L- 
functions that do not include in standard Maple library. In any case we have 
to interval evaluate the expression for zeta or L— function or the closed form of 
value of zeta or L— function. 

4.2 Interval evaluation 
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4.2.1 Riemann zeta 



Consider the series £(s) = ^2 n ^ for complex values of s with Re(s) > 1. 
Interval evaluation of the £(s) can be taken from the result of Backlund (a little 
bit reformulated). The result is based of Euler-Maclaurin summation. 
Proposition (Backlund) Let N be natural > 1. Let s = a + it and let a > 1. 
Let i?2fe be the Bernoulli numbers in even numeration, 

B{N,k,s) = \N 
Then 

where 



n=l 1 s-1 



+ ■■■+ (fffj^s + 1) ■ • • (s + 2fc - 2)iV- Ar - 2 ' £ + 1 . 
C(«)=5(JV-l,s)+B(2V,fc ) s)+ii 2fc , (12) 
2/c - 1 



\R2k 



< 



C7 + 2fc-l 



\B 2 k term of (h 



4.2.2 Dirichlet series 



Let K = Q(\AD) be a real quadratic field with positive integer squarefree D 
and x( n ) — (n) be the Kronecker symbol. Here 

A = 



D, 
AD, 



D = 1 (mod 4) , 
D = 2,3 (mod 4). 



Let A = tt/A, E{xX = C ^-dt, erfc(x) = f /~ e"' 
Proposition [ 

i(i,x) 



tit. 



1 m //i 

^=^ X (n)£;(An 2 ) +^ 

n— 1 n— 1 



er/c(nV^) + i? 



where |i? m | < 



A 3/2 e -J 



4.2.3 Remarks about L~ functions of elliptic curves. 

From the work of A. Wiles, R. Taylor and A. Wiles, and work of F. Diamond it 
is known that (semistable) elliptic curves over E/Q are modular. Knowing the 
modularity of E/Q is equivalent to the existence of a modular form / on Tq(N) 
for some natural value N, which we write / = a n<l n ■ The L— function of E is 
thus given by the Mellin transform of /, L(/, s) = ^a„g n . In particular, the 
behavior of L(E, s) at s = 1 can be deduced from modular properties of E. 
Let -E/Q be a modular elliptic curve and the global minimal model of the E/Q 
has prime conductor I. Let p be a prime and A p (as in (1.4)) be the number 
of points of E p in F p . Then there is exists a modular form / on Tq(1), f = 
J2^Li a n<] n , where a p , p ^ /, equals p + 1 — A p . Under these assumptions it 
seems that results of || gives expressions for L(E, 1) and L (E,s), which can 
evaluate by validated numerics methods. 
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